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Abstract 

The application of optimization techniques derived from the study 
of Euclidean full Steiner Trees to macromolecules like proteins is re¬ 
ported in the present work. We shall use the concept of Euclidean 
Steiner Ratio Function (SRF) as a good Lyapunov function in order 
to perform an elementary stability analysis. 
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1 Introduction 


Nature has followed mathematical principles of structural organization in the 
construction of macromolecular configurations. Our proposal in the present 
work is the modelling of the folded stage of proteins by some combinatorial 
optimization techniques associated to Euclidean full Steiner trees [1]. This 
means that henceforth we take the 3-dimensional Euclidean space E 3 as our 
metric manifold. The analysis to be undertaken can be summarized by the 
trial of obtaining the potential energy minimization of a protein structure 
through the problem of length minimization of an Euclidean Steiner Tree 
[2, 3]. Our fundamental pattern of input points will be given by sets of 
evenly spaced points along a right circular helix of unit radius. We have, 

Pj = (cos ju, sin joj, aju)\ 0<j<n—1, (1.1) 
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where uj is the angular coordinate and 2ira stands for the pitch of the helix. 

We also use the result of Steiner points belonging to another helix of the 
same pitch and smaller radius or 

Sk = (t{u, cc) cos kco, r(u, a) sin kui, akui)\ 1 < k < n — 2, (1.2) 

where 

r(uj,a) = — . = ; A\ = 1 — 2 cos ui. (1.3) 

VMAi + 1) K 

The function r(u, a) above is easily obtained from the requirement of meet¬ 
ing edges at an angle of 120° on each Steiner point. To be rigorous, we 
should write, 


r(ui, a ) = Max [ 1, — _ ] , (1.4) 

v VMAi + 1 ))' 

where the Max above should be understood as a piecewise choice of the 
Maple® software. 


2 Trees of Helical Point Sets 


We now introduce a generalization of the formulae above by thinking on 
subsequences of input and Steiner points, corresponding to non-consecutive 


points. These subsequences 

are of the form: 



^Pmax 

Pj i Pj+rri') Pj +2m •> • • 


(2.1) 

Ismax 

Ski Sk +7715 Sk-\-2mi • • 

• 5 Sk+lsmax-m 

(2.2) 


where Ip, Is are the number of intervals of skipped points before the present 
point on each subsequence and (m — 1) is the number of skipped points. 


We also have: 

0</p< Ipmax = 


n-j-1 


m 


1 ^ Is — ISmax — 


n — k — 2 


m 


(2.3) 


0 < j < m — 1, 1 < k < m 
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and the square brackets stand for the greatest integer value. 


The sequences corresponding to eqs.(l.l) and (1.2) are of course included 
in the scheme above. They are (-Po)i,n-i and (Si)i in _ 2 , respectively. In the 
general case, we can define new sequences of n and n — 2 points instead those 
given by eqs.(l.l) and (1.2). We shall have respectively, 


771—1 m 


Tm — (^k)m,l krnax ■ (2-4) 


j =0 k =1 


The present development is independent of a specific coordinate represen¬ 
tation of the points. If we now assume helical point sets whose points are 
evenly spaced along right circular helices, we get 

Pj+ipm = (cos (j + lpm)u>, sin (j + l P m)uj, a(j + l P m)uj), (2.5) 

Sk+i s m = (r m (w, a) cos (k + lsm)uj, r m (w, a) sin (k + lsm)uj, a(k + l s m)u). 


( 2 . 6 ) 


The function r m (cn, a) is obtained through the same requirement of meeting 
edges at 120° on each Steiner point. We have, analogously, 



(2.7) 


where 


A m = 1 — 2cos(mu;). 


( 2 . 8 ) 


In figure (1) below we show some sequences of input points for n = 23. 
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Figure 1: (A) The sequence n = 23, m = 1, j = 0; (B) The union of the 
subsequences n = 23, m = 2, j = 0 and n = 23, to = 2, j = 1; (C) The union of 
the subsequences n = 23, m = 3, j = 0; n = 23, to = 3, j = 1 and n = 23, to = 3, 
3= 2. 

From eqs.(2.5) and (2.6) and figure (1), we can write for the length of the 
spanning trees 


m— 1 

lsp(m,u>, a) = v m 2 a 2 uj 2 + + 1 . 

j=o 

(2.9) 


n — J — 1 


m 


+(m—1)\/ a 2 u; 2 + Ai + 1. 


The length of the Steiner Trees is then 


lsT(m,u;,a ) 


(1 - r m (uj,a)) 



n — k — 2 


m 


( 2 . 10 ) 
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+ 

+ 


y/m 2 a. 2 u) 2 + r m (u;, a)(A m + 1) . ^ 


k =1 L 


n — k — 2 


m 


2y / m 2 a 2 w 2 + (1 - r m (w,a)) 2 + r m (w,a)(A m + 1). 


After using some useful relations like 

nn —1 


£ 

j=o 


n- j - 1 


m 


— n — m 


£ 

fc=i L 


n — k — 2 


m 


= n — m — 2 


and taking the limit for n m, we get 

lsp(rn,u,a) = n\/ m 2 a 2 u 2 + + 1 


lST ('Pi, oj, a) = n 1 + maw 


A 


A m + 1 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 

(2.14) 


By following the prescriptions for writing the Steiner Ratio, we can write 
for the Steiner Ratio Function of very large helical point sets with points 
evenly spaced along right circular helices 


p(u>,a) 


min 


(m) 


1 + mauj 


/ A m \ 
\ ^4m + l J 


min (m ) (Vm 2 a 2 w 2 + A m + l) 


(2.15) 


where the min process above should be understood in the sense of a piecewise 
function formed by the functions corresponding to the values m = 1,2,3,... 


Eq.(2.15) is our proposal for a Steiner Ratio Function (SRF) [4, 5]. It allows 
for an analytic formulation of the search of the Steiner Ratio which is then 
defined as the minimum of the SRF function, eq.(2.15). Actually, there is a 
further restriction to be imposed on function (2.15) in order to characterize 
it as an useful SRF function. This restriction is that we should consider 
only full Steiner Trees, i.e., non-degenerated Steiner trees in which there 
are exactly three edges meeting at each Steiner point. This restriction can 
be imposed on the spanning trees, by requesting that the angle 9 m (i0,a ) 
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between consecutive edges formed with the points Pj + i p11l as vertices should 
be lesser than 120°. We have 


- < cos 6 m (u), a) 



(-^m + l ) 2 


(2.16) 


In figure (2) below we can see the restrictions corresponding to eq.(2.16), 
for m = 1,2,3. The horizontal line is cos0 m = —1/2. 



cos ft 


-1H 


Figure 2: The restriction to Full Steiner Trees. The figure is a section a = an 
(eq.(2.19)) of the surfaces given by eq.(2.16) corresponding to to = 1,2,3. 

The m = 1 spanning tree is the only one which corresponds to Full Steiner 
trees in a large region of the ^-interval convenient for our work. The other 
trees, m = 2,3 correspond to forbidden regions in the same ^-interval. The 
corresponding Steiner trees to be obtained from the positions of the points 
S'k+isrn and Pj+ipm are necessarily degenerate and should not be taken into 
consideration. Thus, the prescription (2.15) for the SRF function turns into 


p(tv,a) 


1 + aui 



Ai+l 


Max( m ) p m (tv,a) (2.17) 


mm, 


in (m ) (Vm 2 a 2 u}' 2 + A m + l) 


where 



(2.18) 
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The function (2.17) has a global minimum in the point 


and 


(ujr, or) = 7r - arccos 


\/30 


3 ’ 9 (it — arccos |) 


p(u) R , a R ) = ^(3V3 + y/7) = 0.78419037337. 


(2.19) 

( 2 . 20 ) 


For a proof see [4], 

The last value corresponds to the famous main conjecture of ref. [1] about 
the value of the Steiner Ratio in 3-dimensional Euclidean Space. It lead us 
also to think that Nature has solved the problem of energy minimization in 
the organization of intramolecular structure by choosing Steiner Trees as an 
intrinsic part of this structure [6]. 


3 The Stability of Steiner Trees Under Elastic Force 
Deformation 


In the following we continue to work in a M 3 manifold with an Euclidean 
distance definition. Let us now introduce a tree as that of figure (3) 



Figure 3: Geometrical scheme for a Steiner Problem with p = 5. 

There are n input points (position vectors rj) and q = (n — 2)/(p—2) Steiner 
points (position vectors S*,). If q is not an integer number, there is not a 
tree with these n, p values [7]. In figure (3) with p = 5, we assume n to be 
a feasible value. The knowledge of the Steiner Problem tell us that this tree 
structure is not stable since its total length can be reduced by decreasing 
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the number p [ 8 ]. The usual stable Steiner problem corresponds to p = 3. 
In this section we shall give another proof of this fact by exploring the 
concept of a Steiner network with physical interaction among its vertices. 
The structure depicted at figure (3) is a representative of the network which 
models the fundamental interactions inside a biomacromolecule. Let us 
consider the interaction of this structure with similar structures. Let the 
resulting interaction forces as applied to input and Steiner points be fj, f g k , 
respectively and let ls krj be the length of an edge between a Steiner point 
and an input point on its neighbourhood. We have the following identities: 

l S k r j = ^ f? ' (rj-fsJ (3.1) 

Uj 

l S k S h+1 = - 7 - f S k ■ (Sfc+i - S fc ) (3.2) 

a S k 

where aj, ag k , j = 1,2,... ,n, k = 1,2,... ,q stand for the modulus of the 
parallel components to the edges of the resulting forces fj, fs k , respectively. 

The total length of the tree above is 

m— 1 2m —3 3m—5 n q —1 

1 = ^2 ^ 2 ^+ ls 3 r, + • • • + ^ ls q r.j + ^^SfcSfe +1 (3.3) 

j =1 j=m j=2m—2 j=n-m-\- 2 k =1 


From eqs.(2.20), (3.1), we can write the total length in the form 


1 - £ r > - (w + E f< 


j =i 
q -1 


kp— 2fc+l 


£ s * ■ E 

k=2 j=(k—l)p—2k-\-4 

- g , f MSg-i) + £ i 

\ i CLj 

\ V 1 j=n-p-\-2 J 


CLl CLj 

j = l J 

f j ! f s fc , (~ f gfc-i) 
a j a S , fc a Sfe-i 


(3.4) 


We now specialize this set of applied forces at the vertices as being collinear 
with the edges joining them, or 

f j = a jfj II ; f Sfe = a s k fs k II (3-5) 
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where the double vertical stroke means “collinear with the edge” and the 
hat over a letter stands as unit vector. 

We now assume that the forces along the edges are Hooke elastic forces 

f j|| = C ( r J - S k) or fj || = fj = ^(3.6) 

II 1 j ~ &k\\ 

4|| = C'(Sfc +1 - S k ) or 4 h = S k = (3.7) 

where C is the elastic constant. 

The assumption of local equilibrium of these forces lead to the conditions: 



p -1 



Y f j + Si = 0 

3 = 1 

(3.8) 

kp—2k+l 

E ^ 

+ S fc -4-L = 0, k = 2,3,.. 

■,9-1 (3-9) 

j=(k—l)p— 2fc+4 


n 

Y f 3 - 4-1 = 0 • 

j=n—p-\-2 

(3.10) 


This is a set of generalized Fermat problems or Steiner Problems [9]. 

For this equilibrium configuration, eq.(3.3) turns into 

n 

l = ^Z r j • /in- ( 3 - u ) 

3 = 1 

The stability of this equilibrium configuration under a variation of the ap¬ 
plied forces can be tested by 

n 

« = £r; . *4 = 0. (3.12) 

3 =1 

We take cartesian coordinates for the M 3 vectors Yj = ( Xi,yi,Zi ), S k = 
(xs k , ys k i zs k ) and we consider the three independent variations 5 XSk , 5 ySk , 
S Zs in the coordinates of the Steiner points. 
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The corresponding variations in the length of the tree are of the form 


jt s^- x j( r j- s k ) 2 + [{r j -S k )-r j ](x j -x Sk ) 

°x Sh t ~ 0 XSk ' 2^ I| r ,- — Sfcll 3 

3 =1 " 3 


0 (3.13) 


and two other analogous expressions for the variations 5 VSk , S Zs . 


From the arbitrariness of these variations we can write, 


n , 

E — 

3 = 1 


S k ) x [(rj - Sfc) x Tj] 
|l r i - S fc || 3 


= 0 . 


(3.14) 


We can also write 


E 

3 = 1 


• S fc ) 2 - 

l r j - Sfc || 3 


= 0 . 


(3.15) 


We now write the position vectors rj, S k for the configuration depicted at 
figure (3). The points can be taken as evenly spaced along right circular 
helices which radii are 1 and R p , respectively. We have, 

rj = (cos(j - l)w, sin(j - l)w, a{j - l)w) (3.16) 

Sfc = (Rp cos ku :, R p sin ku>, akuj). (3-17) 


Rp(u, a) is a function which can be derived from the equilibrium conditions 
in eqs.(3.7)—(3.9). For p = 3 there is only one solution given by 

„ . . otu 

R 3 (u,a) = . = ; Ai = 1 — 2cosw. (3.18) 

y/Ai(A! + 1 ) 


This solution coincides with eq.(1.3). 

For p > 3, another useful solution could be obtained from the equations: 

cos (fj,ri) = cos (rj,S k ) = cos (S k ,S m ) = 

j,l = l,2 ,---,n; k,m = 1,2,-•• ,q. 


Curiously, Nature has chosen this solution for p = 4 to keep sure of partial 
equilibrium of side chains between the Amide plane conformation in proteins 
[ 6 , 10 , 11 ]. 
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For the configuration given by eqs.(3.15)-(3.16), eq.(3.15) can be written as 


^ ( Tjot) — 0 
3 = 1 

where the geometrical object Tj^p can be written in the coordinates of 
eqs.(3.16) and (3.17) as 

Tjkp = 

[cos 2 (j—1— k)uj— 1— a 2 u> 2 (j— l) 2 ]Rp + a 2 uj 2 [2Rp{j—l)k cos (j—l—k ) u—k 2 } 

[1 + R 2 —2Rp cos(j—\—k)uj + a 2 u) 2 (j— 1— k) 2 } 3 / 2 

(3.19) 

To each £;-value, there will be a term j = k + 1 which dominates the sum 
above. However, we cannot have j = k +1 for p > 3. This can be seen from 
the fact for a vertex {k / 1, q) there are (p — 2) nearest external vertices 
rj. The sequence of their consecutive position vectors is 

r (fc-l)p-2fc+4).> r kp-2k+l (3.20) 

and the requirement j = k + 1 corresponds to an integer p-value only for 
P = 3- 

This p = 3 case which is known to correspond to the most stable problem 
[8] has as a possible configuration the figure (4) below 



Figure 4: The stable structure of the p = 3 Steiner Problem. 
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4 Concluding Remarks 


We have stressed on some past publications that there is a self-consistent 
treatment of the intramolecular organization of biomacromolecules in terms 
of Steiner networks. This representation is able at deriving information con¬ 
cerning its stability and evolution. The supporting facts for stability are now 
well-established and the ideas related to the evolution of macromolecules are 
in their way to be developed and accepted as a preliminary theory of molecu¬ 
lar evolution. The missing subject is a full description of geometric chirality 
and in order to unveil some of its properties, we have proposed to study 
the influence of some proposals for chirality measure on the dynamics of 
optimization problems. These are aimed at studying the structures which 
energy is around the assumed energy of the minimum solution and the vari¬ 
ation process of the chiral properties in the neighbourhood of this minimum. 
We think that this research line is worth of serious scientific work and should 
take advantage of the best efforts of very good scientific researchers for some 
years. 
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